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1. Introduction 



Completion of the knowledge of the generalized nuclear force, which includes not only the 
nucleon-nucleon (AW) interaction but also hyperon-nucleon (YN) and hyperon-hyperon (YY) in- 
teractions, brought the deeper understanding of atomic nuclei, structure of neutron stars and super- 
nova explosions. However it is hard to know the properties of the YN and YY interactions because 
the YN and YY scattering data in free-space are scarce. 

The satisfactory theoretical description of the generalized nuclear force has not yet been pre- 
sented because experimental studies of the YN and YY scatterings are extremely limited due to the 
lack of the hyperon beam. 

Recently a method to extract the AW potential through the Nambu-Bethe-Salpeter (NBS) wave 
function from lattice QCD simulations has been proposed in M\. The obtained potential is found 
to have desirable features, such as attractive well at long and medium distances, and the central 
repulsive core at short distance [Q, The method has been applied to the FAf systems such as AN 
and EN systems [Q, ||, [7|]. Further investigation of baryon-baryon (B-B) interaction have been done 
in the flavor SU (3) symmetric world, and the flavor-spin dependence of S-wave B-B interactions 
are revealed in Ref. [§]. 

In this work, we focus on the S = —2, 1 = B-B system to seek the AA interaction and to 
see the flavor SU (3) breaking effects of the B-B interaction from Lattice QCD simulation. The 
interesting point of this system is that the flavor singlet state is involved. The flavor singlet state is 
known to generate the strong attraction and free from the Pauli blocking at the short range region 
from the view of quark degrees of freedom [p], ||]. For this reason we would like to know how its 
effects appear in the real world. 



2. Formulation 



The (S,I) = (—2,0) B-B state consists of the AA, NE and EE components in terms of low- 
lying baryons. Mass differences of these components are quite small, and it causes a mixture of 
NBS wave function from excited states. In sucn situation the source operator should be optimized 
to extract the energy eigen states through the variational method H [RJ. 

Using diagonalized source operators, we obtain the NBS wave function of AA, NE and ££ 
components for each energy. The transition potential matrix of 3-states coupled channel equation 
can be acquired 



(V 



V'j 



AA 



AA V 
AA 



AA 



V AA V 



/AA 



AA V NE V 



V N hz 

EE / 



(2.1) 



This is an application of an extended HAL QCD method formulated in Ref. [11]. 



2.1 Diagonalization of source operators 

We prepare three wall-source operators at t = to as 

;AA( t \ _ D S \\ \ ^^ t \ rfNSf, \ _ D s 



<f M (to) =P" a p[A a Ap}(t ), jr N *(t ) = P" ap [N a Ep](t ) and ^(t ) = P" ap [L a L p ](t ) (2.2) 
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where the square bracket stands for the 1 = combination of two baryons. The P s is the projection 
operator of spin singlet state. We use the same definition for the interpolating field operators as 
given in Ref. [§]. 

In order to create the energy eigen states, we make a linear combination of three source oper- 
ators as 

>«(0 =v^f AA {t)+^ Z f NE (t)+v^f^(t). (2.3) 
Then we define the effective mass of the operator J? , 



m(At) = — — In 
At 



{j?{t)J{Q)) 



(2.4) 



where At is usually taken to be 1 . The stationary condition of the effective mass against the coeffi- 
cients v leads 

C{t + At)y = e-' n ^- At C{t)y (2.5) 

with the correlation matrix defined as 

^ J {t) = {f\t)J J {Q)) (I,J = AA,NE,IE). (2.6) 

The optimized source operators J?\, J^2 and which strongly couples to the ground, 1st excited 
and 2nd excited state respectively are obtained by solving the generalized eigen-problem. 

2.2 Coupled channel Schrodinger equarion 

The equal-time NBS wave function \j/ B[Bl (r,E) for an energy eigen state with E is extracted 
from the four point function, 

W BiB2 {t-t ,-r') = '£(0\B l {t,x + r)B 2 {t,x)j? E (t ) | 0) , (2.7) 

X 

where J'e is diagonalized wall-source operator which is expected to separate the energy eigen state 
from the other one. 

General form of the Schrodinger equation for the full wave function \jf(f,E) with energy E is 
given in the form with the free Hamiltonian Ho as 

[H -E] ¥ (r,E) = Jd 3 7U(r,7)W,E) (2.8) 

where U (r, r 1 ) is an energy independent but non-local potential. At low energies, it can be expanded 
by the local velocity v = V/ju as U (f, r 1 ) = V(r , V) 8 (r - V) = (V L0 + V NW + V NNL0 + ■ ■ ■ ) 8 (f - ?), 
where N"LO term is of 0(v n ). It is shown in Ref. Jl2| ] that the LO potential Vif) defined above is 
a good approximation of U (r,r') at low kinetic energies between T ~ MeV and T ~ 45 MeV for 
the AW system in quenched QCD. 

Taking the leading order of velocity expansion of non-local interaction kernel and projection 
operators into the AA, NE and ££ components, we have the coupled channel effective Schrodinger 
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Table 1: Lattice parameters and hadron masses in unit of [MeV] are listed. 







Lattice parameters 




p 




csw lattice size 


a [fm] L [fm] 


1.83 


0.13710 


1.7610 16 3 x 32 


0.121 1.93 





Nconf 




«% 


m K 


m N 


m A 


m L 


m= 


Set 1 


700 


0.13760 


875(1) 


916(1) 


1806(3) 


1835(3) 


1841(3) 


1867(2) 


Set 2 


800 


0.13800 


749(1) 


828(1) 


1616(3) 


1671(2) 


1685(2) 


1734(2) 


Set 3 


800 


0.13825 


661(1) 


768(1) 


1482(3) 


1557(3) 


1576(3) 


1640(3) 



equation as 



2}Xnz 

V 2 ft? 
+ 



¥ AA (r,E i )=Y,V AA r (r) ¥ V(r,E i ) 



W ^{r,E i )=Yy m y{r)wy{r,E i ) 
y 



(1 = 0,1,2) 

(y = AA,A^E,EE) 



where asymptotic momenta are related to energies of state as 



Ei = 2Jmi + pf 



rat 



+ (lf + y>n 2 z +q} = 2Jml + kj 



Solving the eq. ( |2.9| ) for the transition potentials, we have 



(2.9) 



(2.10) 



V A % s (r) 



2jUaa 



^(r,^) V^^o) ^ E (r,£o) \ 1 / (V 2 + pi) V^frEo) 



W^Ex) Y N ^(r,Ei) W^{r,E { 



r^(7,E2) W Na (r,E 2 ) Y hh (r,E 2 



v 2 +;? 2j^AA ( 

V {v 2 +pI) x v aa < 



(2.11) 



The other transition potentials are obtained in a similar manner. 



3. Numerical setup 



In this calculation we employ 2 + 1 -flavor full QCD gauge configurations of Japan Lattice 
Data Grid(JLDG)/International Lattice Data Grid(ILDG) [Q. They are generated by the CP- 
PACS and JLQCD Collaborations with a renormalization-group improved gauge action and a non- 
perturbatively 0(a) improved clover quark action at j3 = 6/g 2 = 1.83, corresponding to lattice 
spacings of a = 0.1209 fm jl4|]. We choose three ensembles of the L 3 x T = 16 3 x 32 lattice 
which means the spatial volume of about (2.0 fm) 3 . These ensembles are called as Set 1, 2, and 
3 corresponding to the hopping parameter for light quarks k u j = 0.13760,0.13800 and 0.13825, 
respectively, keeping K s = 0.13710 fixed for the s-quark in all ensemble. Quark propagators are 
calculated from the spatial wall source at to with the Dirichlet boundary condition in temporal 
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Figure 1: Eigen vectors of configuration Set 1. Three panels shows the eigen vector of ground, 1st excited 
and 2nd excited state from the left, respectively. Solid lines stand for the SU (3) Clebsh-Gordan coefficients. 
The "TD" in figures stands for the At in the text. 



direction at t — to = 16. The wall source is placed at 16 different time slices on each of different 
gauge configuration ensembles, in order to enhance the signals, together with the average over 
forward and backward propagations in time. The average over discrete rotations of the cubic group 
is taken for the sink operator, in order to obtain the relative S-wave in the B-B wave function. The 
numerical computation is carried out at KEK supercomputer system, Blue Gene/L. The hadron 
masses are shown in Table [j]. 



4. Result 

Figure [I] shows the eigen vectors of configuration Set 1 . From this figure it is difficult to adopt 
the values in the deep time slice region due to the large noise of the wall-to-wall correlation matrix 
though they are well determined in the shallow time slice region. In this calculation we adopt a set 
of eigen vectors at the shallow time slice but we expect that they could be extrapolate to the deep 
time slice region around t = 9 or 10 because the eigen vectors are not change so much against t. 

The potential matrix V 1 j in Set 1 is shown in Figure fj using the NBS wave function at t — to = 
9. We can see the flavor dependence of the height of repulsive core at short distance region. The ££ 
potential has the strongest repulsive core of these three channels. In off-diagonal parts of potential 
matrix, we find that the strengths of V^zi. and V N ~zz are similar but it of V^nz is much weaker 
than the others. 

In order to compare the results of potential matrix calculated in three configuration sets, we 
transform the potentials from the particle basis to the SU (3) irreducible representation (IR) basis 
as 



V IR = UVU'= V 8 i V 8 V % 2i (4.1) 




where U is an unitary transformation matrix whose explicit form is given in Appendix B in Ref. [j|]. 
The potential matrix in the IR basis is convenient and a good measure of the SU (3) breaking effect 
by comparing three configuration sets because it should be diagonal in the SU (3) symmetric limit. 

In Figure || we compare the results of the potential matrix in the IR basis calculated in different 
configuration sets. We also plot the results of the potential matrix in the SU (3) symmetric limit 
on top of them. Diagonal parts of the potential matrix can be compared with the SU (3) symmetric 
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Figure 2: The potential matrix in Set 1 using the NBS wave function at t = 9 is plotted. The vertical axis is 
for the potential strength in unit of [MeV], and the horizontal for the relative distance between two baryons 
in unit of [fm] . 




Figure 3: Transition potentials in the SU (3) IR basis. Red, blue and green symbol respectively correspond 
to the result of Setl, Set2 and Set3. The results of the SU (3) symmetric configuration with K uc [ s = 0.13710 
are also plotted with brown symbol [||]. 

result. We found the growth of repulstive core in the V 21 potential with decreasing the light quark 
mass. The V1-27 and V%-2j transition potentials are consistent with zero within error bar. On the 
other hand, it is noteworthy that the Vi-g transition potential which is not allowed in the St/ (3) 
symmetric world is appearing. 

5. Conclusion 

We have investigated the (S,I) = (—2,0) B-B state, which is known as the AA, NE and ££ 
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coupled state, from lattice QCD. In order to deal with this complicated system, we have devised 



an extention of HAL QCD method flllQ of construct the interaction potential to a coupled channel 
system. First, we have diagonalized source operators which enable us to obtain the NBS wave 
functions with each energy. Second, we have put all component of the NBS wave function in the 
coupled channel equation for the considering system. Then we could obtain the desired potential 
matrix. Although there have been the challanging problem to determine the eigen vector of source 
operators, we have confirmed that this technique is working well. 

We have searched the SU (3) breaking effects of the (S,I) = (—2,0) BB state by changing the 
light-flavored quark mass. We have found a small transition potential in V\-$ in terms of the SU (3) 
IR basis. Such transition can not be allowed in the SU (3) symmetric world. 

This method could greatly assist us to complete the knowledge of not only the generalized 
nuclear force but also the interaction of hadrons including mesons, baryons and quarks. 
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